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The electron momentum distribution function in the t — J model is studied in 
the framework of slave particle approach. Within the decoupling scheme used in the 
gauge field and related theories, we treat formally phase and amplitude fluctuations 
as well as constraints without further approximations. Our result indicates that the 
electron Fermi surface observed in the high-resolution angle-resolved photoemission 
and inverse photoemission experiments cannot be explained within this framework, 
and the sum rule for the physical electron is not obeyed. A correct scaling behavior 
of the electron momentum distribution function near k ~ kp and k ~ ?>kp in one 
dimension can be reproduced by considering the nonlocal string fields [Z. Y. Weng 
et al., Phys. Rev. B45, 7850 (1992)], but the overall momentum distribution is still 
not correct, at least at the mean field level. 



PACS numbers: 71.45. -d, 75.10. Jm 



I. INTRODUCTION 



The t — J model is one of the simplest models containing the essence of strong correlations 
and its implications for oxide superconductivity still remain an outstanding problem. 
The t — J model was originally introduced as an effective Hamiltonian of the Hubbard 
model in the strong coupling regime, where the on-site Coulomb repulsion U is very large as 
compared with the electron hopping energy t, and therefore the electrons become strongly 
correlated to avoid the double occupancy. In this case, the electron's Hilbert space is severely 
restricted due to this constraint J2a C\ a Ci<T < 1- Anderson [|J and later Zhang and Rice 
have argued strongly that the basic physics of oxide superconductors can be described by 
the t — J model. 

The normal state properties of oxide superconductors exhibit a number of anomalous 
properties in the sense that they do not fit in the conventional Fermi liquid theory [|3],|J. 
Some properties can be interpreted only in terms of a doped Mott insulator . A central 
question in the theory of these strongly correlated systems concerns the nature of the elec- 
tron Fermi surface(EFS) |4j]. The high-resolution angle-resolved photoemission and inverse 
photoemission experiments || demonstrate the existence of a large EFS, with an area con- 
sistent with the band structure calculations. Since the band theory is consistent with the 
Luttinger theorem ||, this means that the EFS area contains 1-5 electron per site, where 
5 is the hole doping concentration. Although the topology of EFS is in general agreement 
with one-electron-band calculations, the Fermi velocity is quite different. This indicates 
that electron correlations renormalize considerably the results obtained in the framework 
of a single-particle treatment. It has recently been shown by small cluster diagonalization 
that in a two-dimensional (2D) square lattice the EFS within the t — J model is consistent 
with Luttinger's theorem |7|]. Monte Carlo simulations for a nearly half- filled 2D Hubbard 
model also support this result ||. Moreover, the electron momentum distribution function 
of a 2D t — J model was studied |§ by using the Luttinger- Jastrow-Gutzwiller variational 



Fermi edge. For one- dimensional (ID) large U limit Hubbard model which is equivalent to 
the t — J model, Ogata and Shiba ]10| obtained the electron momentum distribution func- 



tion by using Lieb-Wu's exact wave function [11]. Their result also shows the existence of 



EFS as well as the singular behavior at k ~ kp and k ~ 3kp in the momentum distribution 



function. Furthermore, Yokoyama and Ogata [12] studied the ID t — J model by using exact 



diagonalization of small systems and found the power- law singularity appearing at kp in the 
momentum distribution function. 

So far, strong correlation effects can be properly taken into account only by numerical 
methods [|7| [TTH , [T^[ ] , such as variational Monte-Carlo technique [JT3], exact cluster diagonaliza- 



tion [jrjj], and various realizations of quantum Monte-Carlo method [[Hifl. Apart from these 
numerical techniques, an analytical approach to the t — J model receiving a great deal of 
attention is the slave particle theory |4jjl6|l, where the electron operator Ci a is presented as 
Cia=a\ficT with a\ as the slave boson and f ia as the fermion, or vice versa. This way the 
non-holonomic constraint 

£<4C iCT <l (!) 

(7 

is converted into a holonomic one 

ata + E/U* = l. ( 2 ) 

(7 

which means a given site cannot be occupied by more than one particle. A new gauge 
degree of freedom must be introduced to incorporate the constraint, which means that the 
slave particle representation should be invariant under a local gauge transformation Oj — > 
a . e iQ(n,t)^ j.^ f ic e t6 ( ri,t \ and all physical quantities should be invariant with respect to this 
transformation. We call such slave particle approach as "conventional". The advantage of 
this formalism is that the charge and spin degrees of freedom of electron may be separated 
at the mean field level, where the elementary charge and spin excitations are called holons 



fluids of spinons and holons represent the same set of electrons, and they must, on average, 
flow together. The decoupling of charge and spin degrees of freedom in the large U limit 
Hubbard model is undoubtedly correct in ID [|Hj], where the charge degrees of freedom of 
the ground state are expressed as a Slater determinant of spinless fermions, while its spin 
degrees of freedom are equivalent to the ID S = \ Heisenberg model, so that charge and 
spin excitations propagate at different velocties. However, the situation is still not clear in 
2D. 

The ID t — J model (the large U limit Hubbard model) behaves like Luttinger liquids 
JHH , |T2"| | . The separation of the spin and charge degrees of freedom is indeed generic to the 
universality class of Luttinger liquids [|TS[]. Anderson [[HJ has hypothesized that the normal 
state of 2D strongly interacting systems relevant for the oxide superconductors should be 
described in terms of the theory of Luttinger liquids. Thus it is interesting to understand 
the global features of the electron momentum distribution function and EFS in the t — J 
model where separation of the spin and charge degrees of freedom might be expected. To our 
knowledge, the global features of the electron momentum distribution function and EFS have 
been obtained only by using numerical techniqus, and have not been studied systematically 
by analytical methods. 

In this paper, we study analytically the electron momentum distribution of the t — J 
model by using the slave particle approach. In Sec. II, we present a formal study of the 
electron momentum distribution function and discuss the sum rule obeyed by it. Under 
the decoupling scheme commonly used in the gauge field JT7J and related theories, we treat 
formally the phase and amplitude fluctuations, constraints and other effects without further 
approximations. The results obtained indicate that the EFS observed in the high-resolution 
angle-resolved photoemission and inverse photoemission experiments f| cannot be explained 
in the conventional slave particle approach. Moreover, the sum rule of the physical electron 
is not obeyed, if the corresponding sum rules are imposed on above particles. Since the 



following Ref. 20, to show that the scaling behavior of the electron momentum distribution 
function near k ~ kp and k ~ 3kp in ID can be reproduced by considering the nonlocal 
string fields in the mean field approximation, but there is still no EFS, and the sum rule for 
the electron number is still violated. If Luttinger theorem is obeyed, the Fermi volume is 
invariant under interactions. Thus our mean field result seems to indicate that the correct 
scaling behavior of the electron momentum distribution does not guarantee the existence of 
EFS (in the sense of global distribution) in the same theoretical framework. Section IV is 
devoted to a summary and discussions on related problems. 



II. FORMAL STUDY OF ELECTRON MOMENTUM DISTRIBUTION 



The slave particle theory can be the slave boson or the slave fermion theories according 



to statistics assigned to spinons and holons. The slave boson formulation [21] is one of 



the popular methods of treating the t — J model. This method, however, does not give 



a good energy of the ground state p2| . For example, in the half-filled case, where the 
t — J model reduces to the antiferromagnetic Heisenberg model, the lowest energy state 
obtained by this method fails to show the expected long-range Neel order, and the energy 
is considerably higher than numerical estimates ||13|| . Also, it does not satisfy the Marshall 
sign rule J23|| . However, this theory provides a spinon Fermi surface []T7| even in the mean 
field approximation. Alternatively, the slave fermion approch naturally gives an ordered 
Neel state at half-filling [p^j , which obeys the Marshall [^3| sign rule. The ground state 
energy obtained in this case is much better than the slave boson case. At the mean field 
level, they are quite different although in princple they should be equivalent to each other. 
The differences may be reduced by going beyond the mean field p5| approximation. In the 



following, we use both approaches to discuss the electron momentum distribution function 
and the sum rule obeyed by it, and a direct comparison of the obtained results is made. 
First consider the slave fermion representation in which the electron operator can be 



with the constraint J2a b\ a bi a + el e « = 1. In this case, the Lagrangian L s f and the partition 
function Z s f of the t — J model in the imaginary time r can be written as 

L sf = £ bld T bu, + £ 4drei + # + £ Ai(et ei + £ b\X« ~ 1), (3) 



H = -t £ j^A- +t E blabipb^bjsia^s - 8 aP 8 jS ) - fi^eUu (4) 

<ij>(T <ij> i 

Z sf = J DeDe ] DbDtfD\e-$ dTL *s {T \ (5) 

where \ is the Lagrangian multiplier on site i, fj, is the chemical potential, and the sum- 
mation < zj > is carried over nearest neighbors. Following the common practice, the boson 
operator b ia keeps track of the spin, while the fermion operator e] keeps track of the charge, 
i.e., they should obey the following sum rules: 

1 c 

5 =< eje. >=— DeD^DbDtfD\eleie-f dTLs f {T \ (6) 
Z 8 f J 

1-5 = J2< b\ a b ia >= -L J2 I DeD^DbDb^D\blb ta e~I dTL ^ T \ (7) 

where 8 is the hole doping concentration, and < ■ ■ ■ > means the thermodynamical average. 
We assume that there is no bose condensation of spinons, which means that the temperature 
24j of the system is T = + . In the gauge field theory |T7| , |26[| , one can introduce the SU(2) 



invariant Hubbard - Stratonovich transformation and decouple the Lagrangian by using the 
auxiliary fields. In the present formal study, there is no need to make any transformation for 
the Lagrangian (3), and we will treat formally phase fluctuations (gauge fields) [[1711 , ampli- 
tude fluctuations ||25|| , constraints et al. For discussing the electron momentum distribution 
function, we define the following Matsubara fermion, boson, and electron Green's functions 

g{Ri -R 3i t-t') = - < T T ei{T)e){T') > 



D a (Ri - R„ T -r') = -< T T b la {r)b) a {r') > 

= --!- / DeDe j DbDb^DXb icT (r)b] a (T')e-I DTL ^ T \ (9) 
Z sf J 

G a (Ri -R J ,r-r') = - < T T C ta {r)C] a {r') > 
= -<T T el(r)e 3 (T')bUr)bl(T')> 
= f DeDe^DbDb j D\el(T)e j (r% (T (T)bUT')e-f dTL °i {T) . (10) 

The spinon-holon scattering contained in eq. (10) is a four particle process, therefore the 
spinon and holon are strongly coupled. At the mean field level, the holons and spinons are 
separated completely. However, in many theoretical frameworks, such as the usual gauge 
theories discussed by many authors |T7|j26[| , where the vertex corrections are ignored, but 



the RPA bubbles are included, the spinons and holons are still strongly coupled by phase 
fluctuations (gauge fields) |T7[], amplitude fluctuations |2Sj and other effects. Nevertheless, 
in all these cases the electron Green's function G a {k, iw n ) can be presented as a convolution 
of the fermion Green's function g(k,iw n ) and boson Green's function D a (k,iw n ) 

G a (k,iw n ) = jj^-^2g(q,iw m )D a (q+ k,iw m + iw n ). (11) 

q M Win 

The coupling of the gauge field to these particles can be strong and a partial resummation 



of digrams has been carried out ||17|| ) but the vertex corrections were neglected, being the 
essence of the decoupling approximation. This is an important but the only approximation 
apart from the sum rules for slave particles (6)-(7) in our formal study. In what follows, 
one will find that many difficulties might appear due to this approximation. The fermion 
and boson Green's functions can be expressed as frequency integrals of fermion and boson 
spectral functions as 

g(k,iw n )= — (12) 



respectively |2"7[] . Substituting eqs. (12) and (13) into eq. (11), we obtain the electron 
Green's function by summing over the Matsubara frequency iw m , 

G a (k,iw n ) = 

I roo dw' f°° dw" . . i\ a t , „,n F (w') +n B (w") 
N g J-oo 2tx J~oo 2it iw n + w' — w" 

where n F (w') and iib{w") are the Fermi and Bose distribution functions, respectively. The 
spectral functions A e (q, w') and A ba (k, w") obey the sum rules coming from the commutation 
relations, 

/CO rjqn 
^A e (q,w) = l, (15) 
-co 2,71 

E/~^M^) = 2. (16) 

The electron's Hilbert space has been severely restricted, but the fermion and boson them- 
selves are not restricted. By an analytic continuation iw n —> w + irj in the electron Green's 
function(14), the electron spectral function A ca (k,w) = —2ImG ca (k,w) can be obtained as 

A ca (k,w) = 

1 f°° dw' 

— J2 —A e (q,w')A ba (q + k,w + w')[n F (w') + n B (w + w')]. (17) 
Therefore, the electron spectral function A ca (k,w) obeys the sum rule 

/°° dm 
^A ca (k,w) = l + S, (18) 
-oo Z7T 

which is less than 2 since an amount of (1 - 5) of the doubly occupied Hilbert space is pushed 
to infinity as U —>■ oo in deriving the t — J model. Thus the spectral function A ca (k, w) only 
describes the lower Hubbard band. In deriving eq. (18), we have used the identities 

/CO f^yj 
—n B {w)A ba (k,w), (19) 
-co 271 



This is because n B (w)A ba (k, w) and n F (w)A e (k, w) can be interpreted as the probability 
functions of state k with energy w for boson and fermion, respectevely. A similar interpretion 
is also valid for the electron spectral function A ca (k,w). Thus the number of electrons in 
state k is obtained by summing over all energies w, weighted by the electron spectral function 

/OO dyj [ ^ 

— n F (w)A ca (k, w) = 1 - 8 - — n ba(k + q)n e (q). (21) 
-oo 2n N qa 

In eq.(21), the first term of the righthand side 1 - 8 is independent of k, and therfore it 
is true for all k states of the entire Brillouin zone. The value of the second term of the 
righthand side is of the order of 8, and hence it is not enough to restore the EFS, i.e., 
the distribution outside the should-be EFS is still of the order 1. Fig. 1 shows the mean 
field electron momentum distribution n c (k) for doping 8 = 0.125 in ID. Beyond the mean 
field approximation, but still within the decoupling scheme (11), there are no important 
corrections for the global features of the electron momentum distribution and EFS, but the 
Fermi velocity will be modified. This is because the essential global features of the electron 
momentum distribution are dominated by the first term of the righthand side, 1 — 8 in eq. 
(21), and the second term of the righthand side in eq. (21), which is of the order of 8, is 
not enough to cancel out 1 — 8 at each k outside the k F state, i.e., k F < k < ir, for k > 0, 
and —tv < k < —k F , for k < 0. Since rib a (k + q) > and n e (q) > 0, and a minus sign 
appears between the first and the second terms of the righthand side in eq. (21), then it is 
impossible to shift the weight 1 — 5 in state k outside the Fermi points into state k' inside 
the Fermi points, i.e., < k' < k F , for k! > 0, and — k F < k' < 0, for k' < 0. The above 
discussions are also true for the 2D case. Therefore, there is no EFS in the standard sense 
within the decoupling scheme (11) in the conventional slave fermion approach. If 8 holes are 
introduced into the half-filled system, one might expect that the total electron number per 
site would be 1 - 8. However, a surprising result is 

k»eW = (l-5) 2 , (22) 



Alternatively, in the slave boson representation, the electron operator can be expressed as 
Cia = fiabl, where f ia is fermion and b\ is slave boson, with the contraint J2& flfia+b\bi = 1. 
In this case, the Lagrangian L s b and the partition function Z s b of the t — J model may be 
written as 

L sb = E flMia + E btdA + H + Hb\h + E flU - 1), (23) 

ia i i & 

H = -t E kblflfja-^flU 

<ij>a ia 

8 a p^s), (24) 

<ij> 

Z sb = J DbDtfDfDfD\e-f drL ° biT \ (25) 

where the fermion operator f} a keeps track of the spin, while the boson operator b\ keeps 
track of the charge, i.e., 

1 c 

S =< b\bi >=-=- DbDtfDfDfD\blbie-f dTL ° b{T \ (26) 
Z s b J 

1 - 8 = E < flU >= ^-E / DbDtfDfDpD\flf la e-f dTL ^\ (27) 

a L sb a J 

We assume that there is no bose condensation of holons. Strictly speaking, this is true above 



the Bose condensation temperature. However, as shown by Lee and Nagaosa [|17| , the in- 
elastic scattering of bosons by the gauge field suppresses significantly the Bose condensation 
temperature. After some formal calculations which are similar to the slave fermion case, we 
obtain 

/OO rfon 
^A ca (k,w) = l + 5, (28) 
-oo ZTC 

n c {k) = 1 -5 + — J2n fa {k + q)n b {q), (29) 

qa 



In comparison with eq. (21), the second term of the righthand side of eq. (29) changes 
sign as an essential difference of the electron momentum distribution function between the 
slave boson and slave fermion representations. However, the value of the second term of the 
righthand side of eq. (29) is also of the order of 5, and it is also not enough to restore the 
EFS. Fig. 2 shows the mean field electron momentum distribution n c (k) for doping 8 = 0.125 
in ID. Beyond the mean field approximation, but still within the decoupling scheme, there 
are no important corrections for the global features of the electron momentum distribution. 
The reason is almost the same as in the conventional slave fermion approach. The essential 
global features of the electron momentum distribution are dominated by the first term of the 
righthand side 1 — 5 in eq. (29). Since nj a (k + q) > and n b (q) > 0, and a plus sign appears 
between the first and the second terms of the righthand side in eq. (29), it is impossible to 
remove those 1 — 5 states beyond the Fermi points. The best situation is that an amount 
of order of 5 is added into each k' state within the Fermi points. The above discussion is 
also valid for the 2D case. Therefore, there is no EFS within the decoupling scheme in the 
conventional slave boson approach. In this case, the sum rule of the electron number is also 
violated as in the slave fermion case. The difference between these two approaches is that 
the electron number is more than the expected value in the slave boson representation, but 
less than it in the slave fermion representation. 

These results indicate that there is no real EFS for the electron momentum distribution 
function n c (k) near k ~ kp within the decoupling scheme in the conventional slave particle 
approach. The sum rule of the physical electron number is violated. The scaling behavior of 
the electron momentum distribution near kp is also not correctly described by this scheme. 
The theory can be applied to both ID and 2D systems. It was proposed |TT| , |2T)| that the low 
energy physics of the t— J model can be described by a theory of fermions and bosons coupled 
by a gauge field. Our results also indicate that the gauge fields (phase fluctuations) discussed 



by many authors [|17|[2t^[28| are not strong enough to restore the EFS under the decoupling 



experiments |J that oxide superconducting materials exhibit an EFS, and the EFS obeys 
Luttinger theorem || , which means the Fermi volume is invariant under interaction. In this 
sense the EFS of strong correlated systems should also be described by an adequate theory 
even in the mean field approximation. In fact, the rearrangements of spin configuations in 



the electron hopping process fLOfl , which are nonlocal effects and beyond the conventional 
slave particle approach, play an essential role to restore EFS for a system of decoupled charge 
and spin degrees of freedom. Thus the electron is not a composite of holon and spinon only, 
and it should include other fields which describe the nonlocal effects. In the next Section, 
we will see that the scaling behavior of n c (k) near k ~ kp and k ~ 3kp can be obtained 
in ID by considering some nonlocal effects. Before going to the next Section, we would 
like to emphasize that the "Fermi Surface" obtained previously |2Pj from the slave boson 
theory is a spinon Fermi surface, not the real EFS. In order to interpret the high-resolution 
angle-resolved photoemission and inverse photoemission experiments ||, the spinon Fermi 
surface is not enough because experiments have shown the EFS for real electrons. 

III. THE SCALING BEHAVIOR OF THE ELECTRON MOMENTUM 

DISTRIBUTION 

FUNCTION AND THE EFFECTS OF NONLOCAL STRING FIELD IN ID 



Interacting ID electron systems generally behave like Luttinger liquids |HJ in which the 
correlation functions have power-law decay with exponents which depend on the interaction 
strength. For the ID Hubbard model, an exact solution was explicitly obtained by Lieb and 
Wu 0]. In the limit of J — > 0, the t — J model is equivalent to the large U limit Hubbard 
model fUJ^JlQH ■ Thus the ID t — J model provides a good test for various approaches. The 
asymptotic forms of some correlation functions as well as the single electron Green's function 
have been obtained by many authors using different approximations |H],|3(]]. In particular, 
Ogata and Shiba |10| obtained the electron momentum distribution function n c (k) by using 



and the correct scaling behavior of n c (k) at k ~ kp and k ~ 3kp. It is well established that 
the Landau Fermi liquid theory breaks down in ID, namely, (1) there is no finite jump of 
the momentum distribution at the Fermi surface; (2) there is no quasi-particle propagation 
and (3) the spin and charge are separated. On the other hand, there is still a well defined 
Fermi surface at kp, as one would expect from the Luttinger theorem. It is remarkable 
that the exact solution of the ID Hubbard model demonstrates explicitly these two aspects 
at the same time. In this sense it is important to check whether both aspects remain in 
any appropriate approximate treatment of the ID model. To our knowledge, such global 
features of the electron momentum distribution function even in ID were obtained only by 
using numerical techniqes. In this Section, we try to study analytically this problem. A CP 1 
boson-soliton approach including the effects of the nonlocal string field to study the large 



U Hubbard model was recently developed by Weng et al. p0| . They have shown that the 
electron is a composite particle of holon, spinon, together with a nonlocal string field. We 
will draw heavily on their results, but we try to make the presentation self-contained. The 
correct scaling form of the electron momentum distribution function can be obtained even 
in the mean field approximation if one considers the nonlocal string field. 

For convenience, we begin with the t — J model, but consider the limit J — > + which is 
a fixed point different from J = 0. In this case, the t — J Hamiltonian may be written as 

H = -t J2 C lCja + h.c. (31) 

<ij>a 

As in the CP 1 boson-soliton approach p0| , the electron operator Ci a can be expressed in 
the slave fermion representation including the effects of nonlocal string fields as 

C l(J = e\e- i ^<^b la G l , (32) 

with the constraint YL C ^leA^ + e l e i — 1- Here rii = ejei, the fermion operator ej keeps track 
of the charge, while the boson operator b\ a keeps track of the spin, e~ l 2 X,i<i n * is the string 



that bi a annihilates the spin a and is a nonlocal phase shift Gi = e l ^^ N ~ Xi ~^- Jl > ini \X i is the 
lattice site), which describes the effects of rearrangements of the spin configurations from site 
i to +00 when one electron was removed or added at site i. As shown by Ogata and Shiba 
TIJ, the "squeezing?' effect, i.e., the rearrangement of spins when holes are squeezed out, is 



crucial in recovering the Fermi surface and getting the correct exponents for the correlation 
functions. The nonlocal "string" proposed in [pOl , will partly take into account this effect. 
In fact the motion of holons is also affected by the rearrangements of the spin configurations 



| 3lfl . We have neglected this second effect since the kinetic energy t is much larger than 
the magnetic energy J. One can check easily that the anticommunication relations are the 
same as for the conventional slave fermion approach. In this case, the Hamiltonian with 
contraints can be written as 



H = -tY^[ib)j> j+ i„e) +1 ej + (-i)b] +l<T b jcr e]e j+1 \ 

~H e] e i + 12 A i ( e i e i + 12 h \° h i° ~ l ) ' ( 33 ) 
j 3 a 

where \i is the chemical potential, and Xj is the Lagrangian multiplier on site j. The factor 
i in the Hamiltonian is very important, which will shift the holon energy spectrum and will 
give rise to the correct scaling form of the electron momentum distribution function. Thus 
we define the new holon operator as 

e j = e i i x >h j , (34) 
and the Hamiltonian may be rewritten as 

H = -tJ2[bl a b i+la hl +1 hi + h.c] - n^2h\hi 

icr i 

+ 5>(/4/^ + 5>iAx-l). (35) 



Following Weng et al. |]20|], we can obtain the asymptotic singular behavior of the electron 



momentum distribution by a similar calculation, and the result is in agreement with theirs. 



factor will contribute an additional power-law decay in the asymptotic single electron Green's 
function, and one may negelect it for simplicity if only interested in the leading contribution. 
The main effect due to the nonlocal string is already present. In the MFA, our situation 
is very similar. Eq. (32) holds exactly in the Neel limit of the spin configuration, and 
therefore the quantum many-body effects are neglected. This does not affect the leading 
long-wavelength behavior. In this paper, however, we are mainly interested in the global 
features of the electron momentum distribution. If Luttinger theorem || is obeyed, the 
Fermi volume is invariant under interaction and a strongly interacting system should also 
show a large EFS even in the mean field approximation. Thus a mean field treatment is a 
useful test for the present approach. 

The mean field approximation to the Hamiltoian(35) amounts to treating A, as a constant, 
independent of position and to decoupling the spinon-holon interaction in a Hartree-like 
approximation by introducing the order paramerters 

X = E < Wk, >, (36) 

<j> =< hfhf +1 >, (37) 

where we have considered two sublattices A and B with i £ A,i + 1 E B. The self-consistent 
equations about X,x,4> , and /U can be obtained by minimizing the free energy. 

Under the same approximation as the one used to obtain the asymptotic single electron 
Green's function(i.e., Gi factor in eq. (32) has been dropped), the single particle electron 
Green's function G ca (k,iw n ) in the present mean field approxmation can be obtained as 

G ca {k,iw n ) = — ^ — 2 - , (38) 

2N ^ lW n + £g_ fc± |(i +( y) - w q 

where 



w k = X - t<fyy k , e k = X- jj,-txik, j k = 2cos(k). 



(39) 



■TT 

A ca (k, w) = — Y,i n B( w q) + n F (e q -k±^{i+s))]S{w + e q - k ±^{i+5) - w q ), (40) 
iV g± 

/°° duo 
—A ca (k,w) = l + 5, (41) 
-oo Z7T 

ra c (fc) = 1 - * - -77 J] n B (w g )n F (£ 9 _ fc±f (42) 
k«#) = (l-5) 2 (43) 

k 

In comparison with eq. (21), we find the holon energy spectrum has been shifted by |(l + 5) 
due to the presence of nonlocal string fields. Fig. 3 shows the electron momentum distribu- 
tion n c (k) in the mean field approximation, including the nonlocal string fields for doping 
5 = 0.125 in ID. The singular behavior of the electron momentum distribution function at 



kp and 3kp is in qualitative agreement with the numerical results p0| , but the sum rule 
for the total electron number is still violated and EFS still does not exist in the present 
framework under the above approximations. Beyond these approximations, the situation is 
not clear yet. If Luttinger theorem [H] is obeyed, then there are no important corrections 
to the global features of the electron momentum distribution function and no EFS beyond 
the mean field approximation in the same theoretical framework. Since the rearrangements 
of quantum spin configurations in the hopping process play an essential role to restore the 
EFS ||10|| , then eq. (32) which holds exactly in the Neel limit of the spin configuration, 



is perhaps not enough to describe all of these quantum spin configurations. Thus our feeling 
is that in this theoretical framework, the correct singular behavior of electron momentum 
distribution may be obtained, but it does not guarantee the existence of EFS. We also note 



that only the asymptotic forms were discussed in the previous works by many authors [30 
while the global features of the momentum distribution function have not been considered 
by their theoretical methods. 



IV. SUMMARY AND DISCUSSIONS 



We have discussed the electron momentum distribution function and the sum rule of 
the electron number in the t — J model by using the slave particle approach. Under the 
decoupling scheme (11) used in the gauge field |l7j and related theories, we have formally 
proved that EFS can not be restored and the sum rule of electron number is violated in 
the framework of conventional slave particle approach. For ID t — J model, the correct 
singularity forms can be reproduced by considering the nonlocal string fields in the special 
mean field approximation, but there is still no EFS and the sum rule is still not satisfied. 

In the present slave particle approach to study t — J model, we have pushed out the upper 
Hubbard band. This means we have neglected the doubly occupied sites in the original slave 
particle approach [^T|: C i<T = a\f i(T + adif-_ a , with a constraint ajdj + J2a fLfia + d\di — 1; 
where a] is the slave boson, f ia is fermion, and di is boson which describes the doubly 
occupied sites, or vice versa. In fact in the original slave particle approach the electron 
spectral function obeys the sum rule of conventional electron, i.e., J2o- I^°oo ^A ca (k, w) = 2, 
but the sum rule of the electron number is still violated. We also find that the sum rule 
of electron number is not violated in the CP 1 representation of electrons, which probably 
means it is a better choice to study the t — J model. 

The decoupling of charge and spin degrees of freedom of electron is undoubtedly correct 
in ID t — J model in the J — ► + limit JlIJ , and ID t — J model behaves as Luttinger liquids 



|10| , |12|| . Anderson |T9[ has hypothesized that the normal state of 2D strongly interacting 
systems relevant for the oxide supercondutors might show generalized Luttinger liquid be- 
havior, and the characteristic separation of charge and spin excitations might be responsible 
for the experimentally observed temperature dependences of the resistivity |32| and Hall 
effect . If these properties can be described in the framework of slave particle approach, 
then we must be out of the difficulties mentioned in Sec. II and Sec. III. A possible way 
to restore the EFS is to try to include the vertex corrections beyond the decoupling scheme 



vertex corrections are neglected, i.e., the decoupling scheme is adopted, the difficulty will 
remain. On the other hand, it is difficulty to introduce vertex corrections in the present form 
of the gauge field theory, because the infrared divergence has not yet been properly handled. 
Another possibility to avoid this difficulty is that one should make a new interpretation of 
the constraint (2) and the physical meaning of spinons and holons. In fact the constraint 
(2) is an operator identity, and one replaces the constraint (2) by eqs. (6) and (7) in the 
slave fermion approach, or eqs. (26) and (27) in the slave boson case. It is not clear that 
this is a correct way of imposing the constraint, i.e., the charge is represented by a fermion, 
while the spin is represented by a boson in the slave fermion approach and vice versa in the 
slave boson version. The crucial point is to implement the local constraint rigorously [J34|| . 
Finally, according to the implications of numerical solutions [ID], the electron is perhaps 



not a composite of holon and spinon only, but it rather should be a composite of holon, 
spinon together with something else. This is because the Bethe-ansatz Lieb-Wu's exact 



wave function []T0[ for ID Hubbard model at large U limit may be written as 

%f(xi,- ■ -,x N ) = (-l) Q det{exp(ik i x Qi )]^(y 1 , ■ ■ -,y M ), (44) 

where the determinant depends only on the coordinates of particles (xq 1 < ■ ■ ■ < xq n ) but 
not on their spins. Thus it is the same as the Slater determinant of spinless fermions with 
momenta kj's. The spin wave function 4>(yi, ■ ■ yu) is the same as the Bethe's exact solution 
of ID Heisenberg spin system. In the second quantization representation, a spinless fermion 
operator may be responsible for the Slater determinant of spinless fermions, but a simple 
fermion or boson operator may be not enough to describe the spin wave function. Perhaps 
it should be a fermion or boson operator together with something else, such as the Jordan- 
Wigner's form [35[] in ID, or a map of the three spin operators of the SU(2) algebra onto a 



couple of canonical boson operators, like the Holstein-Primakoff form ||36|| , are responsible 
for the spin wave function. We believe that a successful theoretical framework must include 



To summarize, we are facing a dilemma: if we would like to split one electron into two 
particles, one keeping track of the charge, the other keeping track of the spin, and impose 
corresponding sum rules on them [eqs. (6)- (7) or eqs. (26)- (27)] and then use the decoupling 
scheme for their expectation values, the resulting electron distribution does not satisfy the 
sum rule and does not show an EFS. In this sense, we have proved in this paper a No- 
Go theorem. The alternative would be to give up the attractively simple interpretation of 
spinons and holons and to look for more complicated charge and spin collective excitations. 

Finally we note that the present slave particle approach is different from those first 
proposed by Barnes P?|] and rediscovered and extended by Coleman Read and Newns 



and Kotliar and Ruckenstein [3D in their works on the mixed valence problem and the 
heavy fermion systems. In their formulation, the spin and charge degrees of freedom are 
not decoupled, where the auxiliary bosons keep track only of the environment by measuring 
the occupation numbers in each of the possible states for electron hopping. Their theory 
describes the properties of a Fermi liquid, and is another story. 
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FIGURES 

FIG. 1. The mean field electron momentum distribution n c (k) for doping 5 = 0.125 in ID in the 
conventional slave fermion approach. The integrated area of the electron momentum distribution 
is 1 -5-5(1 -5) = (1-5) 2 . 

FIG. 2. The mean field electron momentum distribution n c (k) for doping 5 = 0.125 in ID in 
the conventional slave boson approach. The integrated area of the electron momentum distribution 
is 1 - 5 + 5(1 - 5) = 1 - 5 2 

FIG. 3. The mean field electron momentum distribution n c (k) for doping 5 = 0.125 in ID in 
the slave fermion approach including the nonlocal string fields. The integrated area of the electron 
momentum distribution is (1 — 5) 2 . Note that, for 5 = 0.125, 3/cf becomes larger than ir and the 
singularity appears at 3kp-2ir, while -3fcp is less than -it and the singularity appears at -3/cp+27r. 
The values of hp and -3/c^+27r are indicated by arrows. 



